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Abstract 

Previous studies have explored the relationship between termite branch tunnel geometry and 
foraging efficiency in a model simulation in which foraging efficiency, !, for two termite species, 
Coptotermes formosanus Shiraki and Reticulitermes flavipes (Kollar) (Isoptera: Rhinotermitidae), 
was investigated in response to two variables, the probability of tunnel branching {P branch) and the 
probability of tunnel branch termination {Pterm)- It was found that simulated tunnel patterns based 
on empirical data did not have maximum foraging efficiency. We hypothesized that termites 
could increase their foraging efficiency in response to landscape heterogeneity. The present study 
investigated how termites could control the two variables, P branch and Pterm, in response to the 
external environment in terms of tunnel network connectivity. It was found that the best 
simulated strategy for C. formosanus and R. flavipes termites would occur if both P branch and P term 
were increased together. This study provides possible mechanisms for foraging strategies in 
subterranean termites and a baseline for future empirical work. 
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Introduction 



Animals must forage efficiently for food 
resources, or they risk a loss of fitness because 
energy that could be allocated to reproduction 
is used for survival (Werner and Mittelbach, 
1981; Zimmer-Faust, 1987). A wide range of 
strategies could be used to achieve foraging 
efficiency by foragers, depending on the 
environment. Foraging strategies are 
important for understanding species survival 
in terms of foodweb structure and ecosystem 
function and the relationship between species 
biodiversity and ecosystem stability 
(Launchbaugh 1996). Foraging strategies have 
been studied both theoretically (Okubo 1980; 
Shlesinger et al. 1995; Diekmann et al. 2000) 
and experimentally (Cain et al. 1985; Forcardi 
et al. 1996; Gillingham and Bunnell 1989). 

Previous studies have contributed to an 
understanding of foraging behavior (i.e., 
searching for and consuming food) from a 
cost-benefit perspective. However, these 
studies have only investigated foraging in 
open terrestrial environments. The foraging 
strategies of subterranean foragers, such as 
termites, differ in that they are highly 
constrained by the structural design of the 
substrate in which the foragers move (Lee et 
al., 2008a). These constraints cause an 
increase in the number of variables that could 
impact foraging. The foraging behavior of 
subterranean termites has been studied using 
monitoring stations combined with mark- 
release-recapture methods (Grace et al. 1989; 
La Fage et al. 1973; Su et al. 1984). However, 
the information derived from these 
investigations has been insufficient to provide 
an understanding of subterranean termite 
foraging strategy (Gallagher and Jones 2005). 
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Lee et al. (2006, 2007a) studied the foraging 
strategy of subterranean termites and 
constructed a lattice model to simulate termite 
tunnel networks for Coptotermes formosanus 
Shiraki and Reticulitermes flavipes (Kollar). 
The model was based on experimental data 
obtained from homogeneous soil substrates 
without food resources (Su et al. 2004). The 
same model was used to examine the 
relationship between foraging efficiency and 
tunnel network geometry as a way of 
understanding how termites could maximize 
foraging efficiency (Lee et al. 2006, 2007a, 
2007b, 2009a). 

Following these studies, Lee et al. (2009b) 
conducted model simulations using 
experimental values for the probability of 
tunnel branching, Pbmnch, and the probability 
of branch tunnel termination, Pterm- Foraging 
by termites in tunnels that were based on 
experimental values was less efficient than in 
tunnels based on other values of Pbmnch and 
Pterm- These results suggested that termites 
may regulate Pbmnch and P term in order to 
increase the foraging efficiency in response to 
landscape heterogeneity. 

Lee et al. (2008b, c) found evidence for the 
hypothesized role of landscape heterogeneity 
by showing experimentally that termites 
began tunnel excavations only at sites with 
surface irregularities, which are a major factor 
in landscape heterogeneity. 

Because termites have limited energy 
available for digging tunnels and foraging, if 
termites construct too many branching tunnels 
in regions with high heterogeneity, the energy 
available for foraging could be reduced, which 
is likely to lead to a decrease in foraging 
efficiency. Thus an additional constraint 
condition seems to be needed to restrict 
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branching and the occurrence of long 
branched tunnels. This constraint condition 
can be described by the two variables, Pbranch 

and Pterm- 
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1987). Given a network with node set V— {vj, 
V2,..., v n } and a link set between the nodes, its 
Laplacian matrix L is defined as (Merris 
1998): 



The purpose of the present study is to explore 
further the effects of these two variables, 
Pbmnch and P te rm, on foraging efficiency 
because these variables can reveal essential 
features of the termite tunnel network 
(Hedlung and Hederson 1999). Specifically, 
the goal is to understand how Pbranch and P te rm 
can be optimized for higher foraging 
efficiency in relation to tunnel network 
connectivity. 

Methods 

The connectivity of a termite tunnel network 
was used as the constraint condition for values 
of Pbranch and Pterm- The network connectivity 
was characterized by the algebraic 
connectivity, " , defined as the second-smallest 
eigenvalue of the Laplacian matrix of the 
network (Fiedler 1973; Grone and Merris 



degCv,.) J = j 

-1 J* j (adjacent) 

0 , otherwise 



where deg(v z ) denotes the degree of the i-th 
node, which is defined as the number of links 
emerging from the i-th node. The terms / and j 
represent the i-th and j-th nodes, respectively 
(i,j= 1,2,3,..., n). 

The algebraic connectivity, ", reflects how 
well the overall network is connected. In order 
to facilitate an understanding of connectivity, 
three simple networks and their corresponding 
Laplacian matrices are given in Figure 1 as 
examples. The eigenvalues for Figure 1 A are 
listed as {0.0000, 0.4131, 1.1369, 2.3595, 
3.6977, 4.3928}, while the eigenvalues for 




3 


-1 


-1 


-1 


0 


0> 


-1 


3 


0 


-1 


1 


0 


-1 


0 


1 


0 


0 


0 


-1 


-1 


0 


2 


0 


0 


0 


-1 


0 


0 


2 


-1 


0 


0 


0 


0 


-1 


1> 


' 2 


0 


-1 


-1 


0 


0 


0 


: 


0 


-1 


1 


0 


-1 


0 


1 


0 


0 


0 


-1 


-l 


0 


2 


0 


0 


0 


-l 


0 


0 


2 


1 


0 


0 


0 


0 


-1 


1 




-1 


-1 


-1 


0 




-1 


2 


0 


0 


-1 


0 


-1 


0 


1 


0 


0 


0 


-1 


0 


0 


1 


0 


0 


0 


1 


0 


0 


2 


-1 


0 


0 


0 


0 


-1 


1> 



Figure I. Example networks with simple connections and their Laplacian matrices. High quality figures are available online. 
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graphs B and C are {0.0000, 0.2679, 1.0000, 
2.0000, 3.0000, 3.7321} and {0.0000, 0.3249, 
1.0000, 1.4608, 3.0000, 4.2143}, respectively. 
The value of " is highest for network A, 
which is apparent from a visual comparison of 
the three examples. 
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Results 



Connectivity for 1,000 simulated tunnel 
networks of C. formosanus and R. flavipes 
was calculated for different values of P branch 
and Pterm (Figure 3). 



Figure 2 is a schematic representation of a 
tunnel network of C. formosanus. The tunnel 
network was simulated using the lattice model 
proposed in Lee et al. (2009b). The schematic 
consists of 32 nodes, 31 links, and the 
network's Laplacian matrix. 

In the tunnel network, the nodes represent the 
end points of the tunnel segments, defined as 
the line connecting the two closest points of a 
tunnel that did not deviate from the tunnel 
path (Su et al. 2004). In the simulation model, 
when a developing tunnel crosses with other 
tunnels, nodes are not created because 
different tunnels tend to occur at different 
depths below the ground surface. Thus, 
although termite tunnel networks with many 
branching tunnels may appear to have 
complex connectivity among tunnels, their 
actual connectivity can be low. 



For C. formosanus (Figure 3 a), " was higher 
with a higher P term and a lower P branch, while " 
decreased with an increase in Pbranch and a 
decrease in Pterm- Lee et al. (2009b) showed 
that there are three distinct regions in a map of 
foraging efficiency, !, according to Pbranch and 
Pterm- The distinct regions are indicated in 
Figure 3 by the yellow dotted line. In contrast, 
the connectivity map was separated into two 
domains, an upper and lower, as indicated by 
the solid blue line. Silhouette analysis was 
used as the separation method for the " map 
(Rousseeuw 1987). This analysis showed that 
when the " map was partitioned into two 
regions, the optimization score was highest in 
the case of C. formosanus (Table 1). The k- 
means algorithm (Hartigan and Wong 1979) 
was then used to divide the " map with k = 2. 
These divisions are indicated by the blue solid 
lines in Figure 3 (a). The values of " were 
markedly lower in the upper domain than in 
the lower domain. 
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Figure 2. A typical termite tunnel pattern of Coptotermes formosanus, consisting of 32 nodes and 3 1 links and the network's 
Laplacian matrix. High quality figures are available online. 
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Lee et al. (2009b) hypothesized that termites 
control the values of P branch and P te rm in 
response to landscape heterogeneity in order 
to move the red box from region I to region II 
(Figure 3 a). The connectivity map gives a 
constraint condition for moving the red box. 
The best strategy for termites would be to 
increase P branch and P term together in the 
direction of the arrow within the lower 
connectivity domain. 
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Like the " map of C. formosanus, the " map 
for R. flavipes was divided into two domains, 
an upper and lower (Figure 3 b). The red box 
was located at P branch = 0.21 and P term = 0.47. 
According to Lee and Su (2009b), increasing 
Pbranch and decreasing P term could increase 
foraging efficiency. However, in terms of 
network connectivity, increasing both Pbranch 
and P term within the lower " domain in order 
could avoid lower connectivity. 



These data suggest that very low connectivity 
leads to an increase in traveling cost during 
foraging due to numerous tunnel intersections 
that termites have to select depending on their 
direction of movement. On the other hand, 
very high connectivity results in a decrease in 
!, because the short and infrequently 
branching tunnels result in the termites to 
covering an area that it is insufficient to meet 
their foraging needs (see Figures 1 and 2 in 
Lee and Su 2009b). 



Discussion 

Lee et al. (2007a, 2009b), showed that 
simulated termite tunnel patterns, based on 
experimental data obtained from homogenous 
sand substrates, maximize their foraging 
efficiency defined as the ratio of energy gain 
for obtained food to loss for transporting food 
for a given time. In the present study the goal 
was to examine how termites construct their 
tunnel networks in heterogeneous landscape 
and determine if this would optimize their 



Table I. Silhouette score values for each partition 



Species 


Number of Partitions 


2 


3 


4 


5 


6 


7 


8 


9 


C. formosanus 


60.4 


55.8 


55.8 


54.6 


53.8 


56.4 


57.2 


55.4 


R. flavipes 


61.6 


55.4 


54.4 


54.9 


53.5 


53.5 


53.8 


53.3 




0.4, 



031 



(b) 



c 
.a 
a.' 



2021 



I) Higher y 
region 



(II) Lower y 
region 



01 



04 05 06 
P term 



01 

0.1 



02 03 



0.4 0.5 06 
P term 



07 




Figure 3. The map of algebraic connectivity, ", for the branching probability, Pbranch, and the branching termination 
probability, Pterm, for (a) Coptotermes formosanus and (b) Reticulitermes flavipes. Darker shades of gray correspond to lower 
values of connectivity. The solid blue line separates domains of lower and higher " values, and the yellow dotted line divides 
the map into regions with low, middle, and high foraging efficiency, !, respectively. The red box indicates the value of ! for 
empirical tunnel patterns. High quality figures are available online. 
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survival probability in a cost-benefit 
perspective. From the fact that termites begin 
their tunneling behavior only at surface 
irregularity, high heterogeneity is likely to 
cause many tunnel branching because the 
degree of surface irregularity increase with the 
heterogeneity (Lee et al., 2008c). When tunnel 
branching occurs, tunnel intersections are 
generated at the branching site. Termites 
frequently encounter tunnel intersections 
while foraging or moving within their 
networks. The direction chosen by termites at 
tunnel intersections is likely to affect foraging 
efficiency because the path length between 
food resources and the nest can vary 
significantly. Thus, the number of tunnel 
intersections is related to tunnel network 
connectivity, and connectivity is likely to be 
an important factor in foraging efficiency. 

The network patterns were characterized by 
two variables, the probability of tunnel 
branching, P branch, and the probability of 
tunnel branch termination, P te rm-> because these 
variables capture essential features of the 
termite tunnel network (Hedlung and 
Hederson 1999). The simulation showed that 
the two termite species could increase P branch 
and Pterm together to achieve intermediate 
connectivity and thereby improve foraging 
efficiency. 

The results may be inconsistent with field 
observations because there may be other 
constraint conditions such as soil hydrology or 
soil particle size associated with the physical 
environment. Soil characteristics may, in turn, 
interact with resource abundance to affect 
tunnel search patterns. For instance, termites 
could decrease tunnel growth in areas with 
poor soil conditions and increase the extent of 
tunnels in other areas. Although the 
simulation may not provide exact predictions 
of termite behavior in the field, the results of 
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this study provide insights into the foraging 
strategy that could be used by termites to 
improve foraging efficiency. The results also 
suggest directions for future empirical 
investigations of termite foraging strategy. 
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